TPJU 16/2000 
hcp-ph/0012293 



Nuclear effects on the spin-dependent 
structure functions 

ALEKSANDER SOBCZYK0 

Institute of Physics, Jagiellonian University, ul. Reymonta 4, 30-059 Cracow, Poland 

Jerzy Szwee@! 

Institute of Computer Science, Jagiellonian University, ul. Reymonta 4, 30-059 Cracow, Poland 



Abstract 

We address the question how the spin-dependent nucleon structure function 
gi(x,Q 2 ) gets modified when the nucleon is bound inside a nucleus. We analyze 
the influence of nuclear interactions using the A — n model, known to describe well 
the unpolarized effect, and the free polarized parton distributions. The results for 
the neutron in 3 He and proton in 3 H, 7 Li and 19 F are presented, showing significant 
changes in the parton spin distributions and in their moments. Scattering processes 
off polarized Li are suggested which could justify these theoretical calculations and 
shed more light on both nuclear spin structure and short distance QCD. 
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1. The influence of nuclear effects on the nucleon structure functions received enormous 
interest after the measurement by the European Muon Collaboration pQ of the ratio of 
structure functions: 

F N/Fe (r O 2 ) 

' . (") 

F 2 (x, Q 2 ) 

followed by a series of experiments [2j confirming the nontrivial changes of the parton den- 
sities due to the nuclear environment. Many theoretical models have described correctly 
the Bjorken x and A-dependence |3] starting sometimes from quite different assumptions. 
It is therefore not easy to tell what is the underlying dynamics of the "EMC effect" . In 
the search for new tests a polarized version of the effect was proposed some time ago |3]. 
In this paper we take into account theoretical and experimental progress in the knowledge 
of the nucleon spin. 

We will study corrections to the spin-dependent structure function g^(x,Q 2 ) (where 
N denotes proton or neutron) by calculating deviation from unity of the ratio: 



R?(x,Q*) = 9 \ ( ix $ . (1.2) 



Such quantity has become important due to the use of light nuclei (D, 3 He) in the extrac- 
tion of the neutron structure function g»(x, Q 2 ) (SMC 0, E143 6J, E154 [7j, HERMES 
[S]). The definition of g^ A (x, Q 2 ) requires some attention. It is the structure function of 
a single, polarized nucleon inside the nucleus. One should keep in mind that polarizing 
the nucleus in general results in a complicated spin combinations of nucleons and it is 
thus not easy to extract g^ A (x, Q 2 ). Only in selected cases and to certain approximation 
the polarization of the nucleus is equivalent to the polarization of the nucleon in the same 
direction. The nuclei we have chosen below ( 3 He, 3 H and 7 Li, 19 F) are good examples of 
such situation. 

2. The model which describes correctly the x and A-dependence of the unpolarized EMC 
effect was proposed in Ref. Here we extend it to the polarized case. 

We recall conventional picture of nucleus as a system of nucleons and non-nucleonic 
objects: A isobars 10] and excess pions [TT] and assume that polarized deep-inelastic 
scattering may occur either from constituents of the nucleons or from constituents of 
the A isobars. Excess pions are spinless and do not contribute directly in polarized 
scattering. The momentum distributions of non-nucleonic objects come from standard 
nuclear physics calculations f2JE] and are in this sense independent of the model. Such 
construction of the model has an important advantage: we do not have to worry what 
the proposed mechanism does to low energy nuclear physics. 

Another important assumption is that nucleons, A isobars and pions contribute inco- 
herently to the structure function of the nucleus. Thanks to it we can write the "effective 
nucleon" structure function in the nucleus as a sum of convolutions of isolated hadron 
structure functions with momentum distributions taken from nuclear physics 

g^ A (x, Q 2 ) = I* -f N (z)g? (*, Q 2 ) + [* ^f A (z)g? (* Q 2 ) , (2.1) 

Jx Z \Z J Jx Z \Z / 
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where Q is the negative momentum transfer and 

z = A^, a = N,A,n (2.2) 

denotes light-cone {k + = k$ + momentum fraction per nucleon of the interacting nu- 
cleon, A isobar or pion. The distribution functions f a (z) satisfy following sum rules: 

A dzf{z) = l-(n A ), (2.3) 
A dzf A (z) = (n A ), (2.4) 
A dzr(z) = (n*), (2.5) 



o 









as well as the momentum conservation law 

rA 



E/ d* */"(*) = £CO = 1- (2.6) 

OL 

Eqs. (J2.3H2.4jl represent baryon number conservation law, with (wa) and (n n ) defined as 
average numbers of A isobars and excess pions in nucleus, respectively. 
Let us precisely discuss each contribution entering Eq. (|2.1jl . 

2.1. In the parton model our fundamental quantity of interest, the spin-dependent nu- 
cleon structure function g± , can be expressed in terms of differences between the number 
densities of quarks with spin parallel and anti-parallel to the longitudinally polarized 
parent nucleon: 



i N f 

9?W = ^J24{ A( l N (x) + Aq N (x)}, (2.7) 



where 

JV7~\ „N N \^N/\ —N JV 



Aq» (x) = gf - gf , Af (x) = - gf . (2.8) 

In the leading-order (LO) QCD above functions become Q 2 -dependent. In the next-to- 
leading order (NLO) g± can be written as: 



9i 0> Q 2 ) = l -i:e 2 q \Aq N {x,Q 2 )+Aq N {x^) + 



z q 



+ «,(Q 2 ) 



2tt 



AC q ® (Ag^ + Ag w ) + -^AC S ® A<? 



(x,Q 2 ) , (2.9) 



where Ag is polarized gluon density (Ag = g^ — g{), a s (Q 2 ) is QCD running coupling 
constant and Nf denotes the number of active flavors. The spin-dependent Wilson co- 
efficients AC q and AC g in the MS scheme can be found in Ref. [T3j. Convolution ® is 
defined as usual by 

{f®g){x) = £±f( z )g(Zy (2.10) 



In 3 He with the wave function entirely in S state, two protons have opposite spins 
and all the spin is carried by neutron. But in realistic 3 He nucleus higher partial waves 
(S' and D) in the ground state wave function lead to the spin depolarization. Effective 
polarization of neutron in 3 He is estimated in to decrease to V n = (86 ±2)%, whereas 
effective polarization of single proton is V p = (—2.8 ± 0.4)%. Thanks to these quantities 
one can write the spin-dependent structure function g±(x,Q 2 ) from Eq. (|2.1|) as: 

g?(x, Q 2 ) = V n g?(x, Q 2 ) + 2V p g p 1 (x, Q 2 ) (2.11) 

for neutron in 3 He and as: 

g?(x, Q 2 ) = V n g{{x, Q 2 ) + ZP p g?(x, Q 2 ) (2.12) 

for proton in 3 H. In the case of 7 Li we assumed that it consists of 4 He, where two protons 
and two neutrons have opposite spins, and of 3 H, in which we treat the spin depolarization 
like in Eq. (f2~T2l. 

2.2. For nuclei with A > 6 we use an approximate formula jTTj reproducing the effect of 
Fermi motion: 




(2.13) 



kp kp , 

for — < z - 1 < — , 2.14 

win rriN 

and f N (z) = otherwise. The possible corrections to this distribution turn out to be of 
little importance in our problem. The Fermi momenta kp for various nuclei can be found 
in Ref. ^Bj- Here is the nucleon mass. 

The distribution f N {z) for 3 He has been extracted from Ref. |T3] and assumed to 
describe also 3 H. It is worth to note that one can calculate f N {z) from nucleon momen- 
tum distribution pjv(fc) (easily accessible in terms of conventional nuclear theory) using 
approximate relationship |17j : 




(2.15) 



2.3. Even though the pions are spinless and do not directly enter Eq. (|2.1jl . their in- 
fluence comes trough the sum rule (|2.fi|) since baryons share the momentum with pions. 
Effectively, this requires replacement (in Eq. ()2.13|) and next) of z — 1 by z - 1 + (z^), 
where (z n ) is average momentum carried by pions (all distributions are then peaked at 

2.4. The A isobar structure function g^(x, Q 2 ), required in Eq. (|2.1)1 . is not known from 
experiment. A phenomeno logical construction has been presented in Ref. fOj and can 
be extended to the spin-dependent case in a straightforward way. We start from writing 
the valence part of proton and neutron polarized structure functions as: 

gl(x, Q 2 ) = l - (~Au v (x, Q 2 ) + ±Ad v (x, Q 2 )) = ~ (Ja (z, Q 2 ) + -A^x, Q 2 )) (2.16) 



4 



9l(x, Q 2 ) = \ (l^ v (x, Q 2 ) + ^Ad v (x, Q 2 )) = ~ (±A (x, Q 2 ) + ^(x, Q 2 )) (2.17) 

where Ai, I = 0, 1 denotes the contribution in which the valence quark is struck by the 
virtual photon and the remaining (spectator) valence quarks are in spin and isospin / 
configuration. The Ai(x,Q 2 ) can be expressed by the valence quark distributions in the 
proton: 

A (x, Q 2 ) = Au v {x, Q 2 ) - -Ad v {x, Q 2 ), (2.18) 

A 1 (x,Q 2 ) = ^Ad v (x,Q 2 ). (2.19) 

The A structure function is constructed assuming the universality of the functions 
Aj(x, Q 2 ) in ground state baryons. Noting that the spectator valence quarks in the A 
isobar are always in spin-isospin 1 state, one writes: 

(2.20) 

(2.21) 
(2.22) 
(2.23) 

In addition we assume that the sea quarks and gluons remain in the same shape in any of 
the A isobars and neglect Fermi motion effects for the A isobars in our analysis. Hence 
f A (z) has form: 

f A (z) = (n A )5(z-l + (z 7T )). (2.24) 



A++ 


(x,Q 2 ) = 


9 


(x,Q 2 ) 


2 A , 

i = -Ad v 
3 


(x,Q 2 ), 


A+ 

9iv 


(x,Q 2 ) = 


3 1 


[x,Q 2 ) 


= ^AdJ 


x,Q 2 ), 


A , 

9iv 1 


x,Q 2 ) = 


2 -M 

9 u 


x,Q 2 ) 


= l&d v ( 


x,Q 2 ), 


9iv 


(x,Q 2 ) = 


9 


(x,Q 2 ) 


1 A , 

= -Ad v 




(x,Q 2 ). 



3. To demonstrate nuclear effects on the spin-dependent structure functions g{ and g™, 
we have chosen three recent parametrizations of free polarized parton distributions in 
the nucleon at the next-to-leading order (NLO) in the MS scheme: AAC [IB] (called 
AAC-NLO-2), LSS [331 and TBK [20]. For completeness, we also give predictions of 
our model using leading-order (LO) parametrization [18J (called AAC-LO). Nuclear pa- 
rameters required in the calculation are extracted from Ref. [12]. For 3 He and 3 H they 
take the values: (n&) = 0.02, (n^) = 0.05, (z n ) = 0.038, whereas for 7 Li: (tia) = 0.04, 
(n^) = 0.09, (z w ) = 0.069. The Fermi momentum kp for 7 Li is 0.86 fin -1 . In Eq. f!2.1|) 
g A stands for averaged over isospin spin-dependent A isobar structure function. 

The results are presented at Ql = 1 GeV 2 . The first moment of the spin-dependent 
structure function g±(x,Q 2 ) (where (3 = p, p/A ,n ,n/A) is defined as: 

r?(Q 2 )= £ dx gP(x,Q 2 ), (3.1) 

and its values for free proton, free neutron and for proton and neutron in various nuclei 
are presented in all considered parametrizations in Tabled 
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Table 1: First moments rf (Q 2 = 1 GeV 2 ) 





AAC (NLO) 


AAC (LO) 


LSS (NLO) 


TBK (NLO) 


r? 

r p/ 3 H 
r p/ 7 Li 


0.129 
0.110 
0.101 


0.144 
0.123 
0.113 


0.129 
0.109 
0.103 


0.113 
0.095 
0.090 


rn 
1 1 

F n/ 3 He 
1 1 


-0.054 
-0.057 


-0.067 
-0.069 


-0.050 
-0.053 


-0.062 
-0.064 



We expect our model to work for 0.1 < x < 1, since at smaller x possible shadowing 
effects, not included in our calculation, may contribute significantly [21]. The question 
how the first moment of gi is modified can be answered only partially because of this 
limit of applicability of the model. Defining: 

Tl y (Q 2 )= J^dxg^Q 2 ), (3.2) 

in Table [2] we present quantities analogical to those of Table HI but integrated in the 
region 0.1 < x < 1. 



Table 2: First moments Y{ 01 (Q 2 = 1 GeV 2 ) 



AAC (NLO) AAC (LO) LSS (NLO) TBK (NLO) 
0.081 0.097 0.099 0.091 

0.067 0.080 0.082 0.075 

0.063 0.076 0.078 0.071 

r™ 01 -0.024 -0.022 -0.021 -0.020 

r" 7 o" e -0.026 -0.025 -0.024 -0.023 



1 1, o.i 
F P/ 3 H 
1 1, 0.1 
r p/ 7 Li 
1 1. 0.1 



4. We start the discussion from 3 He nucleus, since it is an ideal target to extract the 
neutron structure function g™(x, Q 2 ). We plotted the spin-dependent structure functions 
of neutron in 3 He and free neutron in three different NLO parametrizations in Fig. ^ 
From Table ^ and 121 one reads: 

^^{ 0.96 >x>1.0 

pn/3He \ 0.89 0.1 > X > 1.0 ' 1 1 

what means 11% decrease of the first moment of the spin-dependent structure function 
of neutron in 3 He due to nuclear effects in the range where our models gives predictions. 
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One should remember that in the experimental analyses (like [7J |Hj) a correction for 
the spin depolarization described by Eq. (j2.11|) is often included. That is why it is 
interesting to see what are the corrections predicted by our model not only to the free 
neutron structure function Q 2 ), but also to the function V n gi(x, Q 2 ) + 2V p gi(x, Q 2 ). 
Therefore in Fig. El we presented two ratios: 

Pl nc (x O 2 ) = — HG ( X '^ 2 ) U 2) 



and 



RT\x,Q 2 ) = 91 n( [X n l> (1.3) 
9i{x,Q 2 



in various NLO parametrizations. The ratio Pi He (x, Q 2 ) measures the influence only of 
the A isobars, excess pions and Fermi motion on the spin-dependent structure function of 
neutron in 3 He. If these effects were not important, P| He (x, Q 2 ) would be equal to unity. 
The ratio i?^ He (x, Q 2 ) describes all nuclear effects, including the spin depolarization in 
3 He. We do not present P^ Ue (x,Q 2 ) and R^ Ue (x, Q 2 ) in the whole x region because the 
neutron structure function crosses zero. 

The corrections to the proton structure function g\ (x, Q 2 ) are studied in 3 H and Li 
(we also mention results for 19 F). The 7 Li nucleus is our best example not only for the 
effect which is very pronounced, but also because this nucleus seems to be a realistic 
polarized target. The 3 H is calculated to check the Bjorken sum rule [22J for system with 
A = 3. The spin-dependent structure functions for free proton, proton in 3 H and proton 
in 7 Li are shown in Fig. El in NLO AAC parametrization. We do not plot them in the 
remaining parametrizations, since there is a very small parametrization-dependence in 
our predictions for g v ^ A . In Fig. 0] ratios: 

9i{x,Q 2 ) 

and 

R^Q 2 ) = K]\V (4-5) 

are shown in both NLO and LO AAC parametrization. As compared to 3 He richer content 
of non-nucleonic objects makes the nuclear effect deeper for 7 Li. The resulting corrections 
to the first moments presented in Table 1 and 2 are also considerably larger: 

H^{°- 79 for (46) 
r? \ 0.78 0.1 > x > 1 1 ' 

Having calculated the nuclear effects for both proton in 3 H and neutron in 3 He we 
are able to check the Bjorken sum rule for system with A = 3. The numbers presented 
in Table ^ and 121 show that: 

rf H - vf He „ f Q.91 o>x>i 

Tl-Tl ~ \ 0.88 0.1 >x>l' 1 ' 
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what means 12% reduction in the region where our model is applicable. 

The results for proton in 19 F are very similar to those of 7 Li due to the saturation of 
nuclear parameters (tia), (n n ) and (z n ). This nucleus, less realistic as a polarized target, 
is interesting because of possible application in the hunt for neutralino as a dark matter 
candidate [2~3] . 

One can certainly improve the model presented above. It would be interesting to 
include shadowing effects and extend the model to low x region (x < 0.1). Another 
improvement would be the inclusion of interference terms resulting from N— A transitions. 
Although it is possible to construct the N — A spin dependent structure function in 
analogy to the A, we are unable to extract from the nuclear matter calculation, we are 
basing on ^2], the density f N ~~ A (z) in the nucleus. Approach which attributes all nuclear 
effects (except depolarization Eqs. (|2.11fJ2~T2"J) to the interference term (23] is conceptually 
different from ours: the absolute normalisation of this term is there a free parameter. 

To summarize, we recall the idea how the nucleon spin-dependent structure functions 
get modified due to nuclear environment. The model we have used to present the effect 
has not been chosen by accident. Among other advantages it can be extended from 
unpolarized to polarized version without new, fundamental assumptions. Whereas the 
case of 3 He serves rather only as a warning what size of corrections should one expect when 
extracting the neutron structure function from polarized 3 He target experiments (our 
predictions for neutron in 3 He are within experimental error bars), the 7 Li nucleus seems 
to be more promising. With present experimental techniques one may seriously think of 
deep inelastic polarized lepton - polarized 7 Li scattering or polarized hadron - polarized 
7 Li scattering with direct photon or muon pair production. In all cases the modification 
due to nuclear effects should be measurable. The expected results are interesting for 
both nuclear structure and QCD studies. One should keep in mind that what is usually 
measured in deep inelastic scattering is the asymmetry: 



where Fi is the unpolarized structure function. Since the nuclear effect is similar on both 
the numerator and denominator, one may be misled by the fact that the asymmetry itself 
shows no significant change as compared to the free nucleon case. 
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Figure 1: The NLO neutron structure functions xg™ (x) and xg^ Rc (x) at Q\ — 1 
in various parametrizations 
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Figure 2: The ratios P| He (x) and i?| Hc (x) at Qq — 1 GeV 2 in various NLO paramet 
tions 
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